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Abstract. The huge neutrino pulse that occurs during 
the collapse of a massive stellar core, is expected to con- 
tribute to the origination of a number of isotopes both 
of light chemical elements and heavy ones. In particu- 
lar, evaporation of neutrons from helium nuclei excited 
by neutrino-nuclear inelastic collisions, may result in the 
r-process as it was first discussed by Epstein et al. (1988). 
Here we consider mainly the possibility to obtain the 
considerable amount of neutrons owing to the neutrino 
breakup of helium nuclei. It is shown that, in general, the 
heating of stellar matter due to the neutrino scattering off 
electrons and the heat released from the neutrino-helium 
breakup followed by the thermonuclear reactions should 
be taken into account. On the base of kinetic network, 
using all the important reactions up to Z = 8, the main 
features and the time-dependent character of the neutrino- 
driven neutron flux are investigated. 

The time-dependent densities of free neutrons pro- 
duced in helium breakup, Y n (t), were used to calculate 
the r-process nucleosynthesis with another full kinetic net- 
work for ~ 3200 nuclides. It was found that in the case 
of metal-deficient stars, Z < 0.01 Zq, the resulting den- 
sity of free neutrons seems to be high enough to drive 
the r-process efficiently under favorable conditions. But 
it is impossible to obtain a sufficient amount of heavy 
nuclei in neutrino-induced r-process in a helium shell at 
radii R > R CT « 10 9 cm. We speculate that to make the 
neutrino-induced r-process work efficiently in the shell, 
one has to invoke nonstandard presupernova models in 
which helium hopefully is closer to the collapsed core ow- 
ing, for instance, to a large scale mixing or/and rotation 
and magnetic fields. Apart from this exotic possibility, the 
neutrino-induced nucleosynthesis in the helium shell is cer- 
tainly not strong enough to explain the observed solar r- 
process abundances. 



1. Introduction 

About 20 years ago the gravitational collapse of massive 
stellar cores was recognized to be accompanied by a strong 
neutrino pulse, and a new branch of nucleosynthesis was 
born — the neutrino nucleosynthesis. The general scheme 
of the neutrino nucleosynthesis can be described briefly as 
follows (e.g., Nadyozhin, 1991). Electron neutrinos v e and 
antineutrinos v e can interact with different isotopes (A, Z) 
of chemical elements by means of both charged currents: 



y e + {A,Z) -> (A,Z + l)* +e 
v e + {A, Z) -+{A,Z- 1)* + e 4 
and neutral ones: 
v + {A, Z) (A, Z)* + v' , 



(1) 
(2) 

(3) 



Key words: nuclear reactions 
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where v = v^v^v^v^, v T and v r . Since the typical muon 
and tau neutrino and antineutrino energies are sufficiently 
less than the rest mass of muon and tau-lepton, v^, v r 
and v T interact only by means of neutral current. The 
resulting nuclei are generally produced in highly excited 
states indicated by asterisks. The states decay emitting 
basically neutrons, protons, and a-particles. The prod- 
ucts of reactions (l)-(3) interact both with themselves 
and background nuclei to give rise to a number of issues 
of importance for the origin of chemical elements. The 
existence of the Gamow- Teller and Isobaric- Analog reso- 
nances is crucial for the efficiency of the neutrino nucle- 
osynthesis (see Fuller & Meyer, 1995; Aufderheide et al., 
1994; Panov, 1994, and references therein). 

The beginning of this new issue of nucleosynthesis was 
marked by papers of Domogatsky & Nadyozhin (1977, 
1978, 1980), Domogatsky & Imshennik (1982), and Do- 
mogatsky et al. (1978a, b) devoted to the neutrino-induced 
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production of p-nuclei and a number of light isotopes such 
as 7 Li, 9 Bc, and n B. 

Epstein et al. (1988) (EHC hereafter) put forward an 
elegant idea that, under favorable conditions, the inelas- 
tic scattering of the muon and tau neutrinos off helium 
nuclei in the helium shell would be a good source of neu- 
trons necessary to drive the r-process. They showed that 
the favorable conditions were low metallicity (~ 0.01 Zq) 
and not too large radius of the helium shell (~ 7 x 10 8 cm). 
Woosley et al. (1990) undertook a comprehensive study of 
the neutrino nucleosynthesis and demonstrated, in partic- 
ular, that it was hard to reconcile EHC's statement with 
their supernova models (Woosley & Weaver, 1995) — the 
radius of the helium shell was large in their models and 
the burst of nucleosynthesis stimulated by the shock wave, 
crossing over the helium shell in a few seconds after the 
collapse of a stellar core, was strong enough (according 
to their estimate) to reduce the importance of the previ- 
ous work done by neutrinos. The astrophysical community 
gives credit to Weaver & Woosley for their presupernova 
models which are among the best nowadays. However, it 
would be premature to discard EHC's idea for the only 
reason that it docs not fit these models. Although a ma- 
jor breakthrough in the understanding of the mechanism 
of Type II supernovae has been attained during last years, 
one cannot yet judge with full confidence such details of 
presupernova structure as the characteristics of chemical 
stratification just before and immediately after the on- 
set of the neutrino pulse. The main uncertainties seem 
to be associated with the large scale properties of time- 
dependent convection (Bazan & Arnctt, 1994) and with 
yet poorly studied effects of rotation at terminal stages of 
stellar evolution. It is plausible to imagine, for instance, 
the occurrence of aspherical hydrodynamic flows bring- 
ing helium closer to the stellar center than the radius of 
helium shell in spherically symmetric hydrostatic presu- 
pernova models. 

Though a promising site for the r-process nucleosyn- 
thesis responsible for the solar r-process abundances is ex- 
pected to be the neutrino-driven wind blowing for a few 
seconds from a new-born neutron star (Takahashi et al., 
1994; Qian & Woosley, 1996; Hoffman ct al., 1997), the 
possibility to create heavy elements in the helium shell 
still remains to be of interest. First of all this relates 
to a well developed r-process signature observed recently 
in the extremely metal poor stars (Cowan et al., 1996a, 
1996b; Ryan et al., 1996; Sneden et al., 1996) and possibly 
to some special issues requiring the weak neutron fluxes 
such as the origin of some isotopic anomalies in meteorites 
(Clayton 1989). 

The aim of our study in progress is to specify the con- 
ditions (in terms of helium shell radius, temperature, den- 
sity, and metallicity) favorable for the neutrino-induced 
production of both light isotopes and free neutrons. Here 
preliminary results are reported on free neutron yields rel- 



evant to the r-process, which were already discussed briefly 
in Nadyozhin et al. (1996) and Nadyozhin & Panov (1997). 

2. Properties of the neutrino flux and thermonu- 
clear rates 

We use the neutrino light curve calculated by Nadyozhin 
(1978) for the collapse of iron-oxygen star of 2Mq (Fig. 1). 
The total energy emitted by neutrinos and antincutrinos 
of all the flavors is equal to £ v a = 5.3 x 10 53 erg. 

To obtain the neutrino luminosity in erg s _1 and total 
energy emitted by neutrinos at time t one has to multiply 
l v (t) and E v (t) by £ v p. The neutrino light curve has a 
narrow maximum at t = 0.023 sec and a slowly decaying 
(r ~ 22 sec) long tail. The bulk of the available neutron 
star binding energy is radiated for as long time as ~ 20 sec 
rather than for 3 sec adopted by EHC and Woosley et al. 
(1990). 




} 

-1 1 3 5 7 9 11 13 15 17 19 



t sec 

Fig. 1. The normalized neutrino light curve l v (t) = L v c,/£ I/ p 
and integrated energy of the neutrino flux E vto t{t) = f l„dt, 
E vto t(oo) — 1. (£ v p is the total energy emitted by neutrinos of 
all the flavors) . 

We assume that each neutrino flavor carries away the 
same energy equal to | £ u p and has the Fermi-Dirac ther- 
mal spectrum with zero chemical potential. The mean 
individual energies of electron neutrinos and antineutri- 
nos are taken to be equal to 12 MeV (T ve = 3.81 MeV) 
whereas for muon and tau neutrinos and antineutrinos we 
have chosen a moderate value 25 MeV (T„ MT = 7.94 MeV). 
The cross section of excitation of 4 He by muon and tau 
neutrinos and the branching ratios for proton and neu- 
tron emission were calculated by EHC and Woosley et al. 
(1990) for the neutrino temperatures T„ MT = 4 — 12 MeV. 
We approximated their mean cross section (per helium 
nucleus and per neutrino flavor) by a simple equation 

< thc4„ > = 5.28 x 10 - 43 T 2 exp(-29.4/lV) cm 2 , (4) 
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where T V{1T is in MeV. Equation (4) reflects the threshold 
nature of the process and has an accuracy better than 10% 
within all the temperature range involved (4-12 MeV). 

The number of 4 He nuclei destroyed by neutrinos per 
unit time is given by 



Table 1. The list of nuclear reactions 



*He4 = 



(5) 



5.76 xlCT 2 



(Rs = i?/10 8 cm) , 



where q = 2/3 is the fraction of muon and tau neutrinos 
and antincutrinos in the total neutrino flux; L u p is the to- 
tal neutrino-antincutrino luminosity; R is the helium shell 
radius; and < E v > is the mean energy of muon and tau 
neutrinos and antineutrinos. The number of 4 Hc nuclei in 
unit volume nn e 4 is given by nn e 4 =pN t^Yaei, A^being 
Avogadro's number. So, Yhc4 is the number of 4 He nuclei 
per baryon. The emission rates of neutrons and protons 
(and respectively of 3 He and 3 H) are 



Ho3 



-KY, 



Ho4 • 



Y p = Y m = -6 p r H c4 , (6) 



where b n = 0.471 and b p = 0.516 are the branching ratios 
(Woosley et al., 1990). 

To follow the sequence of thermonuclear reactions 
stimulated by the sources of fresh n, p, 3 He, and 3 H, we 
have used a code connecting 26 nuclides from neutrons and 
protons up to 16 O with a net of about 100 most important 
direct and inverse reactions presented in Tables 1 and 2. 
The reaction rates were taken mostly from Caughlan & 
Fowler (1988) and from the updated version of Thiele- 
mann's original compilation (Thielemann et al., 1986). 

The equations of nuclear kinetics were calculated both 
for constant values of density and temperature and for 
time-dependent temperature induced by the heating of 
stellar matter owing to the neutrino-helium interaction 
and neutrino-electron scattering. 

3. Neutrino heating of stellar matter 

There are two sources of heat worth to be accounted for 
in the helium shell. They are the neutrino scattering off 
electrons and the heat released during the neutrino- helium 
breakup and the subsequent thermonuclear reactions. 



1 H(n, 7 ) 2 H + inv. 


7 Li(a, 7 ) n B 


13 C(n, 7 ) 


14 C 


2 H(p, 7 ) 3 He + inv. 

2 H(p,n)2 1 H 

2 H(D,p) 3 H 


7 Li(a,n) 10 B + inv. 
7 Be(p,7) 8 B + inv. 


13 C(a,n 


16 


14 C(p,n) 


14 N 


7 Be(D,p)2 4 He 


+ inv. 


2 H(D,n) 3 He 


7 Be(T,np)2 4 He 


14 C(p, 7 ) 


15 N 


3 H(p, 7 ) 4 He 


7 Be( 3 He,2p)2 4 He 


13 N( P , 7 ; 


14 


3 H(p,n) 3 He + inv. 


7 Be(a,7) 11 C 
7 Be(a,p) 10 B + inv. 


13 N(n, 7 ; 


14 N 


3 H(D,n) 4 He 


13 N(a,p 


16 Q 


3 H(T,2n) 4 He 
3 He(D,p) 4 He 


9 Be(p, 7 ) 10 B 


14 N( P , 7 ; 


15 Q 


9 Be(p,np)2 4 He a 
9 Be(p,D)2 4 He 


14 N(a, 7 


18p 


3 He(T,np) 4 He 
3 He(T,D) 4 He 


15 n( P , 7 ; 


16 


9 Be(a,n) 12 C b 


15 N(p,a 


12 c 


3 He( 3 He,2p) 4 He 
4 He( 3 He,7) 7 Be + inv. 


10 B(a,n) 13 N 


14 0(n, P ) 

15 o(n, P ; 


14 N 


10 B(p, 7 ) 11 C 


15 N 


4 He(T, 7 ) 7 Li + inv. 
4 He(D, 7 ) 6 Li + inv. 
4 He(T,n) 6 Li + inv. 


n B(p,7) 12 C 


15 0(n,a 


12 C 


11 B(p,n) 11 C + inv. 


18 F(n,a) 


15 N 


n B(a,p) 14 C + inv. 


"Fe"(n, 7 ) 


4 He( 3 He,p) 6 Li + inv. 


n B(a,n) 14 N + inv. 


"Fe"(p, 7 ) 


J-ie(za, 7 J U 


11 TD /„ \ 4tT^ 

d{p,2cx) tie 






4 He(an,7) 9 Be + inv. 


11 C(p, 7 ) 12 N + inv. 






6 Li(a,p) 9 Bc + inv. 


11 C(n, 7 ) 12 C 






6 Li(a, 7 ) 10 B 


11 C(n,2Q) 4 He 






6 Li(p, 7 ) 7 Be 
7 Li(p,n) 7 Be + inv. 


11 C(a,p) 14 N + inv. 






12 C(p, 7 ) 13 N + inv. 






7 Li(p,a) 4 He 
7 Li(D,n)2 4 He 


12 C(n, 7 ) 13 C 






12 C(a, 7 ) 16 O c 






7 Li(T,2n)2 4 He 


13 C(p, 7 ) 14 N 






7 Li( 3 He,np)2 4 He 


13 C(p,n) 13 N + inv. 







a The reaction goes through unstable 9 B^ 2 4 He+p (life time 
~lxl(T 16 s). 

b The rate is taken from Wrean et al. (1994). 

c The rate is multiplied by constant 1.7 (Weaver and Woosley, 



1993). 



Table 2. The weak interaction processes 

4 He(^X) 4 He* 
4 He* ^ 3 H +p 



4 He* 



^ 3 He+n 
,e+)n 
s B-> 2 4 He+e+ 

15 15 N +e 4 



and r neutrinos and antineutrinos, 
the cross sections and branching ratios 
from Woosley et al. (1990). 

r = 0.770 s 
T = 0.011 s 
r = 122 s 

The lifetimes r from Tuli (1990) 



3.1. Neutrino- electron scattering 



The differential cross section of the neutrino scattering off 
electrons at rest {E v » kT), is given by (e.g., Okun, 1982; 
Vogel & Engel, 1989) 



da _ 



~d~£ 



(gv + 9a) 2 + (gv - 9a) 2 (l - + 



(2 2 \ £ 

{9a - 9v) -gi 



where = 2.20 x 10~ 45 cm 2 , E v and e are the energy of 
incident neutrino and the recoil kinetic energy of electron, 
respectively (both in terms of m e c 2 ), and 



9v 



(7) g A = 



2 sin 2 6 W + § 
2 sin 2 6 W - \ 
l 

J 

2 



for v e 
for 

for v e 
for v„ 



(8) 
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Here 9\y is the Weinberg angle (sin 2 6w = 0.23). For 
antineutrinos, one must make the substitution gA — * —gA- 
The kinetic energy e depends on the angle of scattered 
neutrino and varies within the limits 



Qb = N A Yj 



< CHe4i > 



(E* - Bn c4 ) i 
<E vi > 



< £ < £ m = 



2El 
TT2E, 



(9) 



Integrating eda/de over the electron kinetic energy in the 
limits given by Eq. (9) and averaging over the Fermi-Dirac 
neutrino spectrum, we obtain 



£ da' 
He, 



^ /F 4 (0) kTi 



(i 



where 



(10) 



a< = + 79a + W9v9a) 1 , 

h i = \ i 5 9v + g\ + &gvgA) i , 

qi and Tj are the fraction in the total neutrino luminos- 
ity and the Fermi-Dirac temperature of type i neutrinos, 
respectively; F^O) and ^(0) are the Fermi-Dirac func- 
tions (F 4 (0)/F 3 (0) = 4.106). Deriving Eq.(10), we have 
neglected small terms of the order of (m e c 2 / E v ) 2 . 

The rate of energy release per unit mass, Q s is given 

by 



in, I _ d<T 



(11) 



Here X is the mass fraction of hydrogen (in the helium 
shell X= 0). 

Assuming = 1/6, we obtain numerically 



(13) 



where (E* — Bn C 4) i is a branching-ratio averaged value 
and the Fermi-Dirac average energy of individual neutri- 
nos <E vi > = 3.147fcT,. 

To estimate the rate of energy liberation in the ther- 
monuclear reactions, one has to find a difference between 
the total binding energies of composition on two succes- 
sive time steps of calculations, on line with the integration 
of the equations of nuclear kinetics. However, we shall use 
a simple approximation here. The point is that the ther- 
monuclear reactions reassemble rapidly the bulk of the 
neutrino-destroyed helium and convert the rest of it into 
heavier nuclei (such as 7 Li, n B, 12 C and others), thereby 
returning the binding energy of 4 He back in the form of 
heat. To take this effect into account, it is sufficient to 
neglect the binding energy -Brc4 in Eq. (13). Such a proce- 
dure slightly underestimates the energy release since the 
energy liberated in processing of 4 Hc into heavier nuclei 
happens to be omitted. Assuming B^tei = in Eq. (13), 
neglecting the electron neutrinos and antineutrinos, speci- 
fying r Hc = 1/4, E* = 36MeV (Woosley et al., 1990), and 
J^Qi — 2/3, and using Eq. (4) for <<7H e 4i>, we obtain 
numerically 

Qb = 2.557 x 10 18 ^ TV cxp(-29.4/7V) = 



5.017 xlO 17 ^ erg g^s" 1 , 

Ha 



(14) 



Qs 



1.267 xlO 16 ^^- 



(0.4027 T ve + 0.17307V - 0.0657) erg g^s" 1 , (12) 



where T ve and TV arc i n MeV and i?s = i?/10 8 cm. Since 
XV w 2T„ e , the contribution of muon and tau neutrino 
and antincutrino scattering to heating is almost equal to 
that of the electron neutrino and antincutrino. 



3. 2. Neutrino-helium breakup and thermonuclear reactions 

The neutrino-helium breakup supplies heat in the form of 
kinetic energy of the resulting fragments p, n, 3 H, and 3 He, 
which is equal to the difference between the excitation 
energy E* and binding energy, -Bhc4> of 4 He with respect 
to 3 H+p and 3 He+n. Thus, the total release of heat per 
reaction is E* - B Hci with B Hei w 20.6 MeV and w 19.8 
MeV for the proton and neutron channels, respectively. 
The general expression for the rate of energy release per 
unit mass due to the neutrino-helium breakup is given by 



where we made also a substitution T vflT = 7.94 MeV. 
3.3. Neutrino heating 

For simplicity we assume that density p and radius R of 
the helium shell remain unchanged during the heating. 
This is certainly true for the initial narrow peak of the 
neutrino luminosity (Fig. 1). However, the time scale of 
the neutrino light curve tail is comparable to the hydro- 
dynamic time scale in the helium shell, ~ i?/w S ound, and 
such an assumption slightly overestimates the resulting 
increase in temperature. To obtain the temporal behavior 
of temperature, one has to integrate over time the ther- 
modynamic equation of energy: 



E(T,p) = E(T ,p) 



{Qs+QB)dt, 



(15) 



where E is specific energy. For E, we use the Equation 
of State as given by Blinnikov et al. (1996) which takes 
into account ideal electron-positron gas with an arbitrary 
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Fig. 2. Temperature versus time as a result of the neutrino 
heating. The initial temperature Tg(0) = 0.2. 

degree of degeneration, ideal ion gas, and blackbody radi- 
ation. 

Figure 2 shows temperature as a function of time for 
different p as calculated with Eq. (15). For initial temper- 
ature T 9 = 0.2, density p = 3xl0 3 gcm~ 3 , and radius 
R = 5 x 10 8 cm of the helium shell, the heating is notice- 
able but not crucial — T g increases up to « 0.25 in 10 sec 
after the onset of the neutrino flux. The greater density, 
the stronger heating — for instance, at p —3 x 10 5 g cm~ 3 , 
temperature increases up to 0.32. Since the heating scales 
also as R~ 2 , it must be taken into account in the calcula- 
tions of the neutrino nucleosynthesis at radii R<5 x 10 8 cm 
(especially for the carbon and silicon shells). 




Log t (sec) 



Fig. 3. The neutron abundance versus time. 

The curves labelled by numbers correspond to the cases: 

1 Neutron yield from the neutrino-helium inelastic 
scattering (all the nuclear reactions being 
turned off). 

2 No "Fe" and no initial 14 N. 

3 No 14 N, the burning rate of Fe is slowed by 
a factor of 30. 

4 No Fe and the CNO fraction of initial 14 N. 

5 Cosmic initial fraction of Fe and CNO fraction 
of 14 N; Fe burns via 56 Fe(n,7) reaction. 

6 A half of CNO fraction of 14 N, the burning rate 
of Fe is slowed by a factor of 30 (dashed curve) . 

7 Same as Case 6 but CNO fraction of 14 N. 

8 The CNO fraction of 14 N and burning rate of 
Fe is slowed by a factor of 100. 

9 Same as 5 but no burning of Fe (Yp e =const). 
10 Same as 6, but no burning of Fe (dotted curve). 



4. Neutrino-induced production of neutrons 

A number of calculations were made for different (in- 
cluding extreme) assumptions of initial amounts of iron 
and 14 N and of the helium shell radius R, as well. The 
temporal behavior of neutron yield Y n is shown in Fig. 3 
for R = 5xl0 8 cm, T 9 = 0.2, p = 3xl0 3 g cm" 3 , 
Z = O.OIZq, and different assumptions about the neu- 
tron poisons 14 N and 56 Fe. The neutron density n D (t) = 
pN A Y n (t) = 1.81xl0 27 r n (t)cm- 3 . 

Even with no initial 14 N and Fe (Case 2), the neutron 
density falls substantially (by about 4 orders of magnitude 
during the first 10 s) as compared to the amount supplied 
by the neutrino-helium inelastic scattering (Case 1) 1 . This 
happens because the bulk of the neutrino-destroyed 4 He 
proves to be reassembled mostly through 1 H(n,7) 2 H, 
3 Hc(n,p) 3 H, 2 H(D,p) 3 H, 3 He(D,p) 4 He, 3 H(D,n) 4 Hc, 

1 According to Eqs. (5,6), the total amount of neutrons pro- 
duced by the neutrino breakup of 4 He is equal to Y n — 
5.76 xl0~ 2 KY nei /Rl « 2.7 x 10~ 4 for R 8 = 5. 



3 He(T,D) 4 He, and 3 H(T,2n) 4 He reactions. Only as little 
as ~1% of the 4 He breakup products is left unbound back 
into 4 He to synthesize some 7 Li [mostly through the reac- 
tion 4 He(T,7) 7 Li] and other less abundant species such as 
n B, 13 C, and 14 C. 

In Case 4, there is no Fe and the initial fraction of 
14 N is equal to the cosmic fraction of CNO-isotopes mul- 
tiplied by a factor of 0.95 as it follows from the hy- 
drogen burning in well developed CNO-cycle (Yni4 = 
0.95X C no/14 = 0.0523Z). The same initial value of Y mi 
was chosen also for Cases 5, 7, 8, and 9. For the versions 
described by Cases 5-9, the initial abundance of Fe was as- 
sumed to be equal to its cosmic value of 1fc56 = Xp e /56 = 
1.27 x 10~ 3 Z (Anders & Grevese, 1989). The difference be- 
tween Cases 5 and 4 is that in Case 5 the initial abundance 
of Fe was taken equal to its cosmic value but Fe was burn- 
ing at the rate 

4%Ez = -pN A <av> Fenj Y n Y Fe . (16) 

Since the number of the neutrino-produced neutrons 
(Case 1) exceeds considerably the initial number of Fe- 
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nuclei, it takes only about 0.05 sec to burn out virtually 
all Fe. That is why one can hardly notice the difference 
between Cases 4 and 5. This is, of course, an unreal as- 
sumption. The products of Fe burning still continue to 
absorb neutrons. To simulate this effect we calculated two 
more versions with the rate of Fe-burning decreased by 
factors 30 and 100 (Cases 7 and 8, respectively). In other 
words, we have changed Eq. (16) by multiplying its right- 
hand side by an additional factor of 1/30 or 1/100. How- 
ever, in the equation for dY n /dt, the rate of the neutron 
absorption by Fe remained unchanged and equal to the 
unmodified right-hand side of Eq. (16). This means that 
Yfc begin to decrease substantially after capturing 30 or 
100 neutrons. The factor 1/30 in the right-hand side of 
Eq. (16) seems to be a good compromise to simulate the 
consumption of neutrons in the r-process. In Case 9, Fe 
does not burn at all. In Cases 4-6, the neutron abundance 
reaches the magnitude of Y n ~ 10~ 8 and keeps it as long 
as about 5 sec. The duration of the plateau in Cases 7-8 
is much shorter than for Cases 4-6. The comparison of 
Case 10 with Case 6 shows how strongly Y n depends on 
the parameters of Fe-burning at t > 1 sec. 

The initial abundance of 14 N, equal to the total 
amount of CNO isotopes, is certainly the upper limit. In 
moderately massive stars (10-15Mq), the CNO-cycle is 
hardly to be efficient enough to convert all 16 O into 14 N. 
Thus, it is reasonable to take a lower initial abundance 
of 14 N. In Case 6, Y n was calculated with Tni4 initial = 
0.5X C no cosmic/14 (which gives Ym4,/Y Fc ~ 20 — the 
value accepted by EHC) and the rate of Fe-burning de- 
creased by a factor of 30. This set for properties of the 
neutron poisons 14 N and Fe will be considered as stan- 
dard in our further preliminary r-process calculations. 

The resulting Y D proved to be strongly sensitive to the 
helium shell radius R (Fig. 4). For instance for R = 10 9 cm, 
Y n attains a maximum Yn max « 10~ 9 which is about an 
order of magnitude less than for R = 5 x 10 8 cm and lasts 
only for ~ 0.2 sec rather than ~ 5 sec as in Case 6. 

The dependence of Y n on temperature turns out to be 
more complicated (Fig. 5). The solid lines correspond to 
constant temperatures T 9 = 0.2 (Case 6), 0.3, 0.5, and 0.8 
whereas the dashed one shows the effect of neutrino heat- 
ing for initial temperature T 9 = 0.2 with Tg(t) given by 
the lower curve in Fig. 2. The neutron yield falls abruptly 
even for a moderate increase in T. However at tempera- 
tures T g > 0.8, it starts to grow at the very beginning of 
the neutrino flux (t < 0.1 sec). This happens owing to a 
considerable rise in the rates of (p,n), (a,n), and (7,n) re- 
actions. The possibility of the neutrino-driven r-process at 
temperatures as high as T 9 ~ 2 seems to deserve further 
and more detailed investigation. 

5. The r-process 

The time-dependent neutron yield Y n (t) discussed above 
can now be used to calculate the r-process. Instead of 
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Fig. 4. Neutron abundance versus time for different radii of 
helium shell. 
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Fig. 5. Neutron abundance versus time for different tempera- 
tures of helium shell. 

rough approximations for Y n (t) imitating Case 6 in our 
previous calculations (Nadyozhin et al., 1996, Versions I 
and II), we used the time-dependent neutron yield exactly 
as shown for Case 6 in Fig. 3 and its analog for R = 10 9 cm 
(the curve Rs = 10 in Fig. 4) to fulfill the r-process calcu- 
lations. Such an approximate approach to connect the two 
codes, through the neutron abundance Y n (t) seems to be 
adequate for our preliminary study. The accuracy of the 
approximation is briefly discussed in Appendix. 

The r-elements yields were estimated with the aid of 
the kinetic model of nucleosynthesis developed earlier by 
Blinnikov & Panov (1996). In this model, variations in the 
number density n(A, Z) of each nuclide were determined 
taking into account reactions involving neutrons and (3- 
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decays : 

dY AiZ /dt = —X jn (A, Z) ■ Y A , Z - X (A, Z) ■ Y AiZ - 
Y n (t)pN A <a nl (A,Z)v> Y A , Z + 
X jn (A+l,Z)-Y A+ltZ + 
Y n (t) P N A <a ni {A-l,Z)v> Y A -i,z + 
K e {A, Z - 1) • Y AtZ -i - Ke{A, Z) ■ Y AZ + 
Xp{A + i,Z -l)-Pi{A + i,Z -1) 



i=0,l,2,3 



Y 



A+i,Z-l , 



(17) 



dYjdt = 



A,Z 



A 7n (A, Z) -Y n (t)pN A <a n7 (A, Z)v> 



J2 i-\ f j(A,Z)-P t (A,Z) 



=1,2,3 



Y A ,z 



where 



Ya,. 



n(A,Z) 
pN A 



Yn(t) 



n n (t) 
pN A 



(18) 



(19) 



In addition to equations given by Blinnikov & Panov 
(1996), Eqs. (17,18) contain terms describing the beta- 
delayed emission of neutrons, Pi (A, Z) being the probabil- 



ity of emission of i neutrons, 



=0,1,2,3 



Pi = 1). However, 



in the present calculations we take into account the beta- 
delayed emission of only one neutron, namely we assume 
P2 = P3 = 0. Equation (18) for dY D /dt is not used here 
since Y n (t) was taken from Case 6 in Fig. 3. 

Moreover, we have also added terms accounting for the 
capture of electron neutrinos by nuclei with reaction rates 
\ ue {A,Z). The relative role of the ^-capture in Eq. (17) 
turns out to be insignificant for the environments un- 
der consideration. The contribution of the i> e -capture is 
even smaller owing to a high threshold for the neutron- 
rich nuclei, and we neglect it. The neutrino- nuclear cross 
sections were calculated earlier (Panov, 1994; cf. Fuller 
& Meyer, 1995; McLaughlin & Fuller, 1995), with both 
Isobar-Analog and Gamow-Teller resonances taken into 
account. The hope to obtain a noticeable smoothing of 
the abundance curve of the r-process elements due to the 
neutrino-nuclear interaction, seems hardly to be realized 
since the cross sections a ve (A, Z) themselves demonstrate 
a clear even-odd effect that is decreasing, however, in mag- 
nitude when A and Z are increasing. The r-process calcu- 
lations, fulfilled here with account for the capture of the 
electron neutrinos by nuclei, manifested only one notice- 
able effect — the acceleration of the r-process, as it was 
shown by Nadyozhin & Panov (1993). However, in case 
of a larger contribution of the neutrino-capture reactions 
to the r-process or under the conditions of lower neutron 
fluxes involved (owing, for instance, to the exhaustion 
of neutrons at the very end of the r-process), neutrinos 
would be able to decrease the even-odd effect to some ex- 
tent (Panov, in preparation). It should be mentioned that 



our calculations do not take into account the evaporation 
of nucleons from highly excited states of daughter nuclei. 
This effect seems to be only of minor importance for the 
bulk of stable nuclei synthesized in the r-process (Meyer 
et al., 1992) except for the production of p-nuclei (Domo- 
gatsky & Nadyozhin, 1977, 1978). 

It is quite appropriate to notice here that the calcu- 
lations of a v for neutrinos from collapsed stellar cores in 
the approximation of the GT-resonance energy shape with 
the delta function (Fuller & Meyer, 1995; McLaughlin & 
Fuller, 1995), are accurate enough only when < E u > « 
-Egtr- However, in case of typically short-living nuclei 
along the r-process path, the GT-resonance energy proves 
to be rather high (Panov, 1994) and one has to deal with 
a detailed shape of GT-strcngth function in more accurate 
calculations. 

The boundary conditions for the range of the nuclei 
involved into consideration, were specified in the follow- 
ing way: Z m j n =26 and Z max =95, with the minimum and 
maximum A for each Z being determined by the lightest 
stable nucleus from one side, and by the neutron-drip line 
for the heaviest isotope from the other side. Although we 
did not take into account the proton-rich unstable nuclei, 
the full number of nuclides and kinetic equations involved 
in our calculations was as large as about 3200 (the strict 
number depended on mass formulae used). The reaction 
rates entering Eqs. (17) differ by tens orders of magni- 
tude. Thus, the system of equations for nuclear kinetics, 
to be dealt with, is a classical example of a stiff system 
of ordinary differential equations. In our calculations, we 
used one of the most effective methods to solve such a 
stiff system of equations — Gear's method (Gear, 1971). 
Our code has options permitting also to make use of the 
implicit Adams method and the method by Brayton et al. 
(1972). The description of the complete package of solver 
routines and its applications to the r-process calculations, 
can be found in Blinnikov & Panov (1996). A similar full 
network is developed by Goriely & Arnould (1996). 

Nuclear reaction rates were taken mainly from Thiele- 
mann et al. (1986), and mass relations from Hilf et al. 
(1976), though we also considered the dependence of r- 
proccss yields on various mass relations and nuclear re- 
action sets. Beta-decay rates were taken from Aleksankin 
et al. (1981). We evaluated the influence of different sets of 
beta-decay rates, either. The comparison with beta-decay 
rates of Staudt et al. (1990) shows, that for the majority of 
nuclei with even Z, the calculated values are rather close 
to each other. Nevertheless, for some isotopic chains with 
odd Z, the deviations may be rather large. However, all 
these deviations may affect significantly only the isotopic 
yields but not the elemental ones. Here we concentrate 
mainly on the efficiency of the neutrino-induced r-process 
under the astrophysical conditions expected for the helium 
shell while our discussion of influence of different nuclear 
data available is not so detailed; in particular, it concerns 
the beta-decay rates, beta-delayed neutron emission prob- 
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abilities, and the rates of the electron neutrino nuclear 
capture. Results of the calculations are shown in Figs. 6 




-15 1 '■ ' 1 

50 100 150 200 250 

A 

Fig. 6. The r-process abundances = ^ Y(A, Z) as a func- 
tion of atomic mass. Version I : R = 5x 10 cm, duration time 
t — 10 sec (solid line) and 35 sec (dotted line). Version II 
(dashed line): R = 10 9 cm, duration time r = 100 sec. Mass 
formulae are from Hilf et al. (1976). 



and 7 for Versions I and II and for various mass formu- 
lae. For Version I, independent of the mass formulae used, 
the nuclei from A = 80 through A = 130 and then up to 
A = 200 were synthesized by the 10-th and 35-th second, 
respectively. Under the conditions of Version II it is pos- 
sible also to form some r-process nuclei but not far from 
the iron peak. 

The differences in the shape of the abundance curves 
based on different data sets are not very high, except for 
the region of low A, where the initial nuclei burn out a 
little bit slower when the data of Hilf et al. (1976) and 
reactions rates of Thielemann (1986) are used instead of 
combining the data of Janecke & Eynon (1976) with re- 
actions rates of Panov (1995, 1997) (see Figs. 6 and 7). 

In our network calculations, we have taken into ac- 
count the emission of beta-delayed neutrons to study 
whether they are able to reduce the amplitude of the odd- 
even effect. The smoothing of abundances due to the emis- 
sion of beta-delayed neutrons proved, however, to be small 
in our case, since the neutron flux is not very high and r- 
process nucleosynthesis occurs not far from the stable nu- 
clei having small values of the beta-delayed neutron emis- 
sion probabilities. This effect depends mostly upon the 
changes in r-process path after the break-down of nuclear 
statistical equilibrium and usage of the kinetic network of 
nucleosynthesis (Howard et al., 1993). 



As it can be seen from Figs. 6 and 7, in the region 
of A > 150 the odd-even effect is a little bit more pro- 
nounced - maybe because of the boundary effects such, for 
example, as a step- like shape of the stable nuclei boundary 
(Panov & Blinnikov, in preparation). 




Fig. 7. Same as in Fig. 6, but Mass formulae are taken from 
Janecke & Eynon (1976), and nuclear reaction rates from 
Panov (1995, 1997). 

An additional smoothing of the odd-even effect could 
occur due to emission of the beta-delayed neutrons dur- 
ing subsequent beta-decays into stable isotopes. But under 
the astrophysical conditions discussed here, the r-process 
path lies not far from beta-stable nuclides and only a frac- 
tion of the nuclei formed in sufficient quantity have signifi- 
cant values of beta-delayed neutron emission probabilities. 
That is why in the present case the smoothing of the odd- 
even effect due to delayed neutrons happens to be small 
and cannot change significantly the calculated abundance 
curves. 

6. Discussion and Conclusions 

In this work we began a systematic study of the conditions 
favorable for the neutrino-induced production of free neu- 
trons in the amounts sufficient to form heavy elements 
through the r-process in stellar shells surrounding the col- 
lapsed core — a source of the intense neutrino flux. The 
problem in question is described by four main parameters: 
radius of the shell R (the neutrino flux scales as R~ 2 ), 
density p, and temperature T in the shell, and chemical 
composition of the shell. The latter parameter splits in 
two important ones — the main constituent (helium, car- 
bon, silicon etc.) and the admixture of the neutron poison 
species like 14 N and 56 Fe, the latter being both a raw ma- 
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terial for the r-process and a neutron poison as well. We 
concentrate here on the possibility of the r-process in a 
helium shell consuming neutrons from the breakup of 4 He 
by neutrinos. 

We explored the dependence of the neutron yield, 
Y n , on R varying within 4xl0 8 - 4xl0 9 cm (Fig. 4) for 
constant values of temperature, density, and mctallicity: 
T 9 = 0.2, p = 3xl0 3 g cm~ 3 , and Z= O.O1Z . The cal- 
culations of the r-process with Y n (t) given by the curves 
for R = 5 x 10 8 and 10 9 cm in Fig. 4 show that there exists 
a critical radius of the helium shell R cr ~ 10 9 cm, above 
which the neutrino flux fails to stimulate the r-process 
efficiently. 

The calculations of Y n for different temperatures of 
the helium shell (Fig. 5) demonstrated a strong sensitivity 
of Yn to temperature. It was also shown that in accu- 
rate calculations one has to take into account the heat- 
ing of helium shell by neutrinos, especially if the density 
is large enough (p > 10 5 g cm -3 ). For the initial temper- 
ature T g > 0.5, Y n begins to increase with temperature, 
first at the very onset of the neutrino flux and then pro- 
gressively at longer times (compare curves T 9 = 0.5 and 
0.8 in Fig. 5). If the initial temperature of the helium shell 
happens to be as large as Tg s=s 2, one can expect a certain 
increase in the critical radius R cr - From one side, such tem- 
peratures prevent the neutrino-destroyed helium to be re- 
assembled and thereby support Y n at a higher level. From 
the other side, at such a high temperature, thermonuclear 
reactions (7,p), (7,n), (p,n), (a,n) become fast enough to 
contribute themselves to Y n even when the neutrino flux 
is absent (Woosley et al., 1990). 

It is worth noticing that the shock wave, passing across 
the helium shell in a few seconds after the collapse of stel- 
lar core, heats matter up to Tg = 2-3 and reprocesses the 
products of the neutrino-induced nucleosynthesis. This ef- 
fect deserves more detailed study. Here we would like only 
to mention an interesting chance (see also Domogatsky et 
al., 1978b) given by a possibility of collapses not leading to 
standard powerful supernova explosions. This chance can 
be high if, e.g., the neutron star formation rate is higher 
than the rate of collapsing supernova outbursts. Unfortu- 
nately, the uncertainties of observations do not permit to 
prove this to be true, yet they do not reject this chance, 
either. It is conceivable, that the rotation of a presuper- 
nova plays the major role during those almost silent col- 
lapses. So, a large mass of the material would be kept at 
the radius lower than i? cr by the angular momentum con- 
servation. Later, as suggested by Bisnovatyi-Kogan (1970, 
1980), the magnetic field, enhanced by the differential ro- 
tation, can eject the outer layers (already irradiated by 
the neutrino flux). It is interesting to note that the lat- 
est computations (Bisnovatyi-Kogan et al., 1995) fail to 
produce an explosion on the supernova energy scale. So, 
a powerful shock wave, able to dissociate the r-process 
nucleosynthesis yields, is not formed, but the mechanism 
seems to be able to eject (with relatively low velocity) 



an amount of r-elements sufficient for enrichment of the 
interstellar matter. 

However, it is not yet clear quantitatively how impor- 
tant the possible changes both in the rates A/3, A„ e , and 
their ratios could be for the neutrino nucleosynthesis is- 
sues. These changes are connected with a possible depen- 
dence of A/3 and \ ve on temperature and density owing 
to, for example, the contribution of excited nuclear states 
and the Pauli principle for electrons. 

The conditions chosen by us, in particular R = 
5 x 10 8 cm, appear to be privileged in a certain sense. Here 
the supernova neutrino pulse leads to a neutron yield al- 
ready sufficient to support the r-process through helium 
breakup, but it is not yet adequate for the direct formation 
of heavy nuclei during the r-process due to an increase in 
the neutron excess by the electron neutrino absorption and 
the neutrino-induced evaporation of neutrons from heavy 
nuclei (Domogatsky & Nadyozhin, 1978). Figure 8 shows 
the influence of the neutrino capture by heavy nuclei on 
the r-process. In order to single out this effect, the emis- 
sion of beta-delayed neutrons was not included in these 
two runs. We see only an acceleration of the r-process but 
it is hard to find any additional smoothing of the odd-even 
effect. 

When R decreases, the rate A„ e increases, and the 
growth of the element synthesis due to neutrino capture 
relative to the beta-decays at the waiting points can lead 
finally to smoothing of the element yield curve. This issue, 
however, needs further investigation. 




9 ' 5 75 85 95 105 115 125 135 

A 

Fig. 8. The influence of the neutrino capture by heavy nuclei 
on the r-process for duration r = 10 sec (the capture is on — 
solid line, off — dashed line; Y n (t) for Case 6). 

Recently, it was shown that the extremely metal poor 
stars are made of matter enriched with heavy elements 
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synthesized in a well developed r-process (Cowan et al., 
1996a, 1996b; Ryan et al., 1996; Sneden et al., 1996). In 
the near future, we are going to extend our calculations 
to metallicities as low as Z= O.OOIZq and even less to 
estimate the resulting increase in R cr . 

In general, we conclude that it is very difficult to obtain 
the neutrino-induced r-process in a helium shell at radii 
R > R cr ~ 10 9 cm, in a good agreement with Woosley 
et al. (1990). The current presupernova models have radii 
of the helium shell at least several times larger than R cr 
specified above. Nevertheless, we continue to search for the 
physical conditions appropriate for the neutrino-driven r- 
process, not binding ourselves with the question of com- 
patibility with the current presupernova models and hav- 
ing in mind probable modifications of the models involved, 
in particular, large scale convection, rotation, and mag- 
netic fields. 
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Appendix: a selfconsistent solution for Y n 

We use two nuclear kinetic codes. One deals with 101 reac- 
tions involving 24 light nuclides from 1 H and n through 
18 F (Code L hereafter), listed in Tables 1 and 2. In ad- 
dition, it takes into account the absorption of neutrons 
and protons by 56 Fe. Another code (Code H), designed 
for studying the r-process, connects more than 3000 heav- 
ier nuclides in a network of neutron-capture reactions 
and beta-decays. The two codes are physically connected 
through the exchange of neutrons (the proton-capture by 
Fe in Code L is virtually unimportant at rather low tem- 
peratures considered here). The r-process converts initial 
Fe into heavier nuclei which continue to absorb neutrons 
at a rate that changes as the r-process nucleosynthesis 
"wave" reaches successively heavier and heavier nuclei. 
Figure 3 shows the results of various assumptions con- 
cerning the rate of the neutron absorption in the r-process, 
made in Code L, to supply Code H with the neutron abun- 
dance Y n (t) that would permit to calculate the r-process 



without appealing to Code L any more. Case 6 was chosen 
as a compromise. 

There are two ways to check the accuracy of such an 
approximation. The most straightforward way would be to 
unify the two networks in one code. This time-consuming 
work has been planned for the future. Another way is to 
iterate Y n (t) function as described below. 

First, one calculates the r-process using an approxi- 
mate solution for time-dependent neutron density Y® (t) 
from Code L instead of Eq. (18). The calculations provide 
the rate of the neutron density change, required by the r- 
process Y®(t) and given by the right-hand side of Eq. (18). 
Then, using this Y® one begins to run Code L with slightly 
modified differential equation for Y n : 

dY n /dt = B + Y n f(t), (20) 

where f°(t) = Y°(t)/Y°(t) and B denotes all the neutron 
absorbing and emitting reactions given in Tables 1 and 
2 except for the reaction Fe(n, 7). In Case 6, the second 
term in the right-hand side of Eq. (20) equals exactly to 
the unmodified right-hand side of Eq. (16) as described in 
section 4. As a result of this Code L run one gets an it- 
erated function Y*(t). Entering Code H again, now with 
Y*(t), one obtains Y^(t) that gives the second approxi- 
mation Y^{t) for Y n (t) after running Code L and so on 

Figure 9 shows the results of such an iterative fairly 
fast-converging procedure. It takes 4-5 iterations to get 
the final selfconsistent solution shown in Fig. 9 with stars 
nonconnected by line. The thin line underlying the stars 
represents the last but one iteration that is almost indis- 
cernible from the final solution. 

From Fig. 9, one can conclude that Case 6 ensures a 
reasonable approximation for Y n (t) as long as t <5 s. How- 
ever, for longer time (t > 10 s), Case 6 overestimates Y D (t) 
by an order of magnitude — the r-process consumes about 
90% of neutrons provided by the neutrino breakup of 4 He. 
We remind that in Case 6 initial Fe proves to be totally 
burnt by t w Is. Therefore for longer times, Case 6 actu- 
ally neglects the absorption of neutrons by the r-process. 
Using the selfconsistent solution for Y n (t) (stars in Fig. 9) 
instead of that for Case 6 (dashed curve in Fig. 9), we have 
recalculated the r-process for times t = 10 and 35 s as 
shown in Fig. 10. Although the r-process yields changed 
at t = 10 s but not so significantly as to alter our con- 
clusions in the main text of the paper. At t = 35 s, the 
changes are crucial — owing to a shortage in neutrons, 
the r-process gets virtually exhausted and fails to form 
the abundance peak at A w 200. To make the r-process as 
far as A = 200 and beyond under physical conditions dis- 
cussed here (T 9 = 0.2,i?=5xl0 8 cm) one needs apply to 
lower metallicities (Z < 0.01 Zq). For the selfconsistent 
solution, the total number of neutrons consumed in the 
r-process per seed 56 Fe nucleus is equal to about 50. 
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Fig. 9. The iterative procedure converting the initial approx- 
imate Y n (t) for Case 6 (dashed line) into a final selfconsistent 
solution (stars). 
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Fig. 10. Same as Version I in Fig. 6 but for selfconsistent Y n (t) 
from Fig. 9. 
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